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1.  INTRODUCTION 


Formation  in  flight  occurs  when  a  radioactive  rare  gas  (contained  within  an  inert  gas  such  as  air), 
after  entering  a  conduit  through  a  high-efficiency  filter,  exhibits  decay  giving  rise  to  a  steady 
production  of  a  certain  number  of  daughter  atoms  per  unit  of  volume.  Unlike  the  radioactive  gas,  the 
daughter  atoms  adhere  to  the  wall,  and  hence,  are  lost  by  diffusion.  This  dispersion  of  atoms  may  be 
considered  to  be  an  atomic  aerosol  (Tan  and  Hsu  1971)  since  it  has  the  same  property  as  submicron 
aerosols  in  that  they  can  be  collected  at  a  surface.  Formation  in  flight  diffusion  equations  have  been 
derived  for  cylindrical  tubes  (Tan  and  Hsu  1971;  Berezhnoi  and  Kirichenko  1964)  and  for  flat 
channels  (Berezhnoi  and  Kirichenko  1964);  all  cases  considered  were  for  slug  (uniform-velocity 
profile)  or  for  Poiseuille  pipe  flow  (parabolic-velocity  profile). 

The  situation  of  a  vortex  or  swirling  flows  involving  just  the  fluid  phase  has  been  investigated  by 
various  authors  (Fromm  1963;  Pao  <967;  Textor  1968;  Munson  1980).  A  well-defined  configuration 
of  a  confined  axially  degenerating  vortex  flow  was  introduced  by  Lavan,  Nielsen,  and  Fejer  (1969). 
More  recently,  Tung  and  Soo  (1973)  treated  the  case  of  swirling  pipe  flow  of  suspensions  for  the  same 
defined  geometry  of  Lavcn,  Nielsen  and  Fejer  (1969). 

The  present  study  applies  the  results  of  the  flow  system  of  Laven,  Nielsen,  and  Fejer  (1969)  and 
Tung  and  Soo  (1973)  (see  Figure  1)  to  the  case  of  a  decaying  radioactive  gas  under  gravitational 
effects;  the  products  of  disintegration  of  the  gas  are  filtered  out  at  far  upstream  (z  =  ><»),  but  are  again 
generated  by  radioactive  decay  of  a  flowing  inert  gas  along  the  entire  length  of  the  cylindrical  tube. 
The  radioelements  diffuse  radially  to  the  tube  walls  where  they  decay  into  other  radioelements. 

The  fluid  phase  is  assumed  incompressible  and  fully  developed  at  both  far  upstream  and 
downstream  positions  from  the  junc.ure  of  the  two  pipes.  A  forced  vortex  is  generated  by  the  rotating 
pipe,  and  it  degenerates  or  dissipates  as  the  flow  passes  through  the  stationary  duct  due  to  frictioa  At 
very  far  downstream,  the  swirl  vanishes  and  the  axial  velocity  reverts  back  to  the  fully  developed 
laminar  parabolic  profile.  The  solutions  of  the  fluid  phase  velocities  can  be  used  to  solve  for  the 
particulate  distribution  governed  by  the  diffusion  equation. 

A  vertical  configuration  will  b r.  assumed  such  that  the  axial  flow  direction  is  coincident  with 
positive  gravity  having  vertical  symmetry. 
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Figure  1.  Coordinate  System  and  Boundary  Values. 


2.  FORMULATION 


The  motion  of  the  fluid  phase  is  given  by  solving  the  continuity  and  the  Navier-Stokes  equations 
for  a  steady,  laminar,  incompressible,  and  axisymmetric  swirling  flow.  The  solution  is  given  in  terms 
of  the  radial,  tangential,  and  axial  components  of  the  fluid  velocity  and  pressure  (Laven,  Nielsen,  and 
Fejer  1969;  Tung  and  Soo  1973)  such  that, 


p  dr 


±  »<»]♦£ 

r  dr  dz2 


9v  MV  3v 


13.  A  d2v 

—  T-(rv)  +^T 
r  dr  dz2 


dw  dw  1  dp 
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~  (ru)  +  JL  (rw)  =  0.  (^) 
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Equations  1  through  3  represent  the  momentum  equations  of  the  fluid  phase  in  the  radial, 
tangential,  and  axial  directions,  respectively.  Equation  4  is  the  continuity  equation  for  the  fluid  phase. 
The  panicle  diffusion  equation  with  gravity  and  generation  can  be  written  as: 


(7) 
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where  ft  is  the  constant  angular  vel  icity  of  the  rotating  pipe.w  is  the  mean  axial  flow  velocity,  and  R0 
is  the  radius  of  the  pipe.  The  flow  configuration  together  with  the  boundary  conditions,  is  show’n  in 
Figure  1 .  The  pipe  section  at  the  left  of  the  r-axis  rotates  with  constant  angular  velocity  ft  and 
generates  the  swirl,  while  the  pipe  section  at  the  right  is  held  stationary.  At  z  =  0,  the  two  pipes  join 
smoothly  with  no  effect  on  the  flow  . 

The  results  for  the  fluid  phase  are  expressed  in  terms  of  nondimensional  variables  by  defining. 


R  =  r/R0,  Z  -  z/R0,.'¥ 
w  =  U  -  JL,  V  -  v/(P  ft) ; 

vv/2  w/2  '  ' 

C  -  taR?!'Vw,p  -p/(pp02ft2);  00) 


where  is  the  stream  function  at  the  wall.  Furthermore,  the  Reynolds  number  is  given  by 


Rey  - 


2wR0 

i 

v 


01) 


4 


and  the  swirl  ratio  is  given  by  the  ratio  of  the  tangential  velocity  of  the  rotating  pipe  to  the  mean  axial 
velocity,  i.e., 
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ORJ 
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(12) 


The  classical  numerical  approach  which  converts  the  velocity  components  into  stream  function  'P, 
vorticity  £,  and  circulation  T  was  chosen  because  of  its  well-established  stability.  Similarity 
transformation  (Kidd  and  Farris  1969)  was  not  used  because  it  is  not  always  possible  to  obtain  a 
physically  realistic  velocity  distribution. 

Numerical  calculations  were  carried  out  for  various  Reynolds  number  Rey  and  swirl  ratio  S 
combinations;  they  were  chosen  to  exhibit  the  various  effects  of  the  swirl  on  the  flow  field.  In  some 
cases,  the  axial  velocity  can  reach  zero  or  even  negative  velocity  (reversed  flow)  at  some  downstream 
axial  positions.  Correctness  of  the  numerical  method  and  computer  program  was  checked  against  the 
results  of  Laven,  Nielsen,  and  Fejer  (1969)  and  Tung  and  Soo  (1973).  The  Reynolds  numbers  and 
swirl  ratios  were  chosen  to  achieve  rapid  convergence  and  consistency  with  Laven,  Nielsen,  and  Fejer 
(1969)  and  Tung  and  Soo  (1973). 

The  diffusion  equation.  Equation  5,  was  also  cast  into  nondimensional  form  by  defining 

C*-Dpc/(QR0 2).  03) 


along  with  the  aid  of  the  diffusive  Peclet  number,  P,  =  R0w/Dp,  and  the  gravity  flow  parameter, 
o  =  gR0/(FDp)  such  that. 
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where  the  mapping  function,  given  by  Laven,  Nielsen,  and  Fejer  (1969),  was  used  to  transform  the 
z-axis  into  the  x-axis.  Also,  F  is  the  inverse  of  relaxation  time  constant  defined  as, 


F  «  9ji/(2 a2c). 


(15) 


where  p  and  c  are  the  viscosity  of  the  fluid  and  material  density  of  the  panicle  phase,  and  a  is  the 
radius  of  an  atomic  panicle. 

Since  it  was  assumed  that  the  fluid  phase  is  not  affected  by  the  presence  of  the  solid  atomic 
panicles  (Soo  1969),  the  fluid  velocities  obtained  from  die  fluid  phase  relations  can  be  directly 
substituted  into  the  paniculate  diffusion  equation.  After  the  particulate  concentrations  were  found  at 
each  respective  axial  and  radial  position  throughout  the  flow  field,  the  bulk  concentration  was 
determined  to  estimate  the  deposition  from  far  upstream  to  downstream  axial  positions. 

Defining  the  F(X)  as  the  ratio  of  the  total  particle  flux  over  a  cross-section  at  distance  X  from  the 
far  upstream  position  (at  z  =  -«>)  to  the  rate  of  formation  of  the  radioelements  in  the  same  element  of 
the  tube,  we  have 


F(X)  - 


1 

X 


r  1  * 

WC  RdR. 

J0 


(16) 


The  relations  for  the  fluid  phase  were  cast  into  explicit  finite  difference  molecules  (Southwell 
1940;  Ames  1972)  and  solved  by  numerical  relaxation  yielding  results  consistent  with  Laven,  Nielsen, 
and  Fejer  (1969)  and  Tung  and  Soo  (1973).  As  a  result  of  the  axial  symmetry  of  the  flow  field,  only 
the  flow  in  the  rectangular  region, 

D=  [(RX)  \0ZRZ1  mdOZXZl],  (17) 

need  be  considered.  On  D,  a  network  of  parallel  mesh  points  (41  x  41)  was  uniformly  spaced  by  an 
amount  of  A  in  both  R  and  X  directions.  With  the  initial  assumed  values,  r*,  £*,  and  'F*  are  solved 
iteratively  by  a  procedure  consisting  of  making  sweeps  of  the  interior  mesh  points  until  a  convergence 
criteria  given  by  Forsythe  and  Wasow  (1967)  is  satisfied. 
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Since  Equation  14  is  parabolic,  and  due  to  the  universally  stable  implicit  method  (Hombeck 
1975),  the  particulate  diffusion  equation  was  cast  into  the  implicit  finite  difference  discretization  form 
as, 


where 
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(19) 


(20) 


(21) 


(22) 


Convergence  was  satisfied  for  all  Peclet  and  gravity  parameters  considered  in  Equation  18  with  respect 
to  the  given  mesh  sizes  considered  in  the  analysis.  The  sparse  and  tri -diagonal  system  resulting  from 
the  finite  difference  expressions  for  the  diffusion  equation  was  solved  by  Gaussian  elimination  with 
full  and  partial  pivoting  aided  by  the  Crout  reduction  technique. 

3.  RESULTS  AND  DISCUSSION 

The  results  of  the  fluid  phase  are  virtually  the  same  as  Laven,  Nielsen,  and  Fejer  (1969)  and 
Tung  and  Soo  (1973),  the  only  difference  being  with  respect  to  precision  which  was  within  0.01%  for 
most  of  the  comparable  calculations. 
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Figure  2  depicts  the  axial  variation  (from  z  =  -°°  equivalent  to  x  =  0  to  1)  of  the  bulk 
concentration  of  the  disintegration  product  formation  of  particles  in  varying  gravity  fields  at  low  and 
moderate  Peclet  numbers.  At  the  same  gravitational  influence,  increasing  the  inertial  forces  over  the 
diffusive  forces  causes  an  increase  in  penetration;  also,  increasing  the  gravity  field  at  a  fixed  Peclet 
number  increases  the  penetration  or  inhibits  wall  deposition. 

As  seen  in  Figure  3,  increasing  the  gravity  field  at  a  fixed  inertial  to  diffusive  force  ratio  caused  a 
smaller  concentration  profile;  whereas,  at  a  fixed  gravity  field,  increasing  the  Peclet  number  diminishes 
the  axial  distribution  of  the  concentration  profile. 

The  effects  of  varying  the  Reynolds  number  and  swirl  ratio  at  a  given  gravity  flow  and  Peclet 
parameter  (viz.,  Rey  e{l,  4,  10,  20},  Swirl  Number  e{  12.  7,  8},  respectively)  did  not  significantly 
alter  the  deposition,  perhaps  due  to  the  restrictive  range  of  the  Reynolds  and  Swirl  parameters. 
However,  because  of  the  need  to  be  consistent  with  Laven,  Nielsen,  and  Fejer  (1969)  and  Tung  Soo 
(1973),  and  for  the  desire  to  limit  computer  time,  the  above  limits  were  set  for  this  investigation. 

The  bulk  concentration  of  particles, 


c  (z)  *  2 71 


J  wcrdr 

0 _ 

W7t  Rg 


is  related  to  F(X)  by: 


/)  n  *  t)  R 

c(z)  -  XF(X)  Z-L  -  XF(X)  Z-l  Pt 


w 


This  shows  clearly  that  decreasing  the  particle  diffusion  coefficient  increases  the  particle  bulk 
concentration  at  each  point  in  the  tube  if  all  other  parameters  are  unchanged  as  exemplified  in  Figures 
2  and  3;  increasing  the  gravity  parameter  has  a  similar  effect. 
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Figure  3.  Axial  Distribution  of  the  Centerline  Concentration  at  Different  Gravity  Fields  (Re y  =  20  and  S 
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c 

D' 

F 

g 

P 

P 

P. 

Q 

r 

R 

Ro 

Rey 

S 

u 

U 

v 

V 

w 

W 

w 

X 

r 

z 

Z 

a 

pl.  P2,  p3 

P 

v 


•  Radius  of  particle 

-  Particle  concentration 

-  Dimensionless  particle  concentration 

-  Material  density  of  particle  phase 

-  Particle  diffusivity 

-  Inverse  of  relaxation  time 

-  Local  acceleration  of  gravity 

-  Dimensionless  pressure 

-  Fluid  static  pressure 

-  Diffusive  Peclet  number 

•  Rate  of  particle  generation  per  unit  of  volume 

-  Radial  coordinate 

-  Dimensionless  radial  coordinate,  r/R, 

-  Pipe  radius 

•  Reynolds  number 

-  Swirl  ratio 

-  Radial  velocity 

-  Dimensionless  radial  velocity,  uVw/2 

-  Tangential  velocity 

•  Dimensionless  tangential  velocity 

-  Axial  velocity 

•  Dimensionless  axial  velocity 

-  Mean  axial  velocity 

-  Transformation  variable  [0,1] 

-  First  derivative  of  X  with  respect  to  Z 

•  Axial  coordinate 

-  Dimensionless  axial  coordinate 

-  Electrostatic  charge  parameter 

-  Finite  difference  coefficients 

-  Viscosity  of  the  fluid 

•  Fluid  kinematic  viscosity 
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A R 
AX 

r 

l 

o 


¥ 

r*.  C,  *f* 


n 


•  Radial  change 

-  Axial  change 

-  Circulation 

-  Vorticity 

-  Gravity  flow  parameter 

-  Stream  function 

-  Dimensionless  circulation,  vorticity,  and  stream  functions,  respectively 

-  Constant  angular  velocity 
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